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Abstract. Daganzo's criticisms of second-order fluid approximations of traffic flow [C. Daganzo, Transpn. 
Res. B. 29, 277-286 (1995)] and Aw and Rascle's proposal how to overcome them [A. Aw and M. Rascle, 
SIAM J. Appl. Math. 60, 916-938 (2000)] have stimulated an intensive scientific activity in the field of 
traffic modeling. Here, we will revisit their arguments and the interpretations behind them. We will start 
by analyzing the linear stability of traffic models, which is a widely established approach to study the 
ability of traffic models to describe emergent traffic jams. Besides deriving a collection of useful formulas 
for stability analyses, the main attention is put on the characteristic speeds, which are related to the group 
velocities of the linearized model equations. Most macroscopic traffic models with a dynamic velocity 
equation appear to predict two characteristic speeds, one of which is faster than the average velocity. 
This has been claimed to constitute a theoretical inconsistency. We will carefully discuss arguments for 
and against this view. In particular, we will shed some new light on the problem by comparing Payne's 
macroscopic traffic model with the Aw-Rascle model and macroscopic with microscopic traffic models. 
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PACS. 89.40.Bb Land transportation 
83.60.Wc Flow instabilities 



45.70.Vn Granular models of complex systems; traffic flow 



1 Introduction 

Understanding traffic congestion has puzzled not only 
traffic engineers, but also a large number of physicists [UEi 
3,4 . Scientists have been particularly interested in emer- 
gent traffic jams, which are related to instabilities in the 
traffic flow. Such instabilities have been found in empirical 
data [5], but also in recent experiments E|. 

The theoretical analysis is usually done by computer 
simulation or by linear stability analysis. Both techniques 
have been used since the early days of traffic engineering 
[7] and traffic physics 8,9 . Here, we will perform the anal- 
ysis for macroscopic and microscopic models in parallel, as 
there should be a correspondence between the properties 
of both kinds of models. In contrast to previous publica- 
tions, the analysis of macroscopic traffic equations is done 
for a model that considers a dependence of the optimal ve- 
locity function and the traffic pressure on the average ve- 
locity, not only the density. Such a dependence results for 
models which represent vehicle interactions realistically, 
taking into account a velocity-dependent safety distance 
[10] . This is, for example, important to avoid accidents, 
and it changes the instability conditions significantly (see 
Sec. El). 

Besides determining the stability threshold, a partic- 
ular focus will be put on the calculation of the group ve- 
locities of the partial differential equations underlying the 
macroscopic traffic model (see Sec. 13. 2j) . For clarity, the 



definition of the group velocities will be compared with 
those of phase velocities and of characteristic speeds. All 
three definitions describe propagation processes of waves. 
It will be shown, that they lead to identical results under 
certain circumstances, but not necessarily so. 



Furthermore, we will derive conditions under which 
one of the group velocities is greater than the average 
velocity. In Sec. [2J we will shortly summarize the main 
points of the controversial discussion that this observa- 
tion has triggered. We will also address Daganzo's other 
criticisms of second-order fluid approximations of traffic 
flow [IT] . After the formal analysis in Sec. El Sec. 2] will 
be dedicated to a careful discussion of the results. In par- 
ticular, we will analyze different conceivable reasons for 
characteristic speeds faster than the vehicle speeds: (1) 
artifacts due to approximations underlying second-order 
macroscopic traffic models, (2) indirect long-range forward 
interactions with followers on a circular road, (3) the def- 
inition of the propagation speed of perturbations, (4) the 
variability of vehicle velocities, (5) the interpretation of 
characteristic speeds. Since characteristic speeds are pri- 
marily perceived as a problem of second-order macroscopic 
traffic models, in Sec. El we will compare them with the 
group velocities predicted by microscopic traffic models. 
Finally, we will summarize our results in Sec. El 
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2 Summary of the Controversy regarding 
Second-Order Traffic Flow Models 

In the area of macroscopic traffic flow modeling, it is com- 
mon to formulate equations for the vehicle density p(x, t) 
as a function of space x and time t and for the aver- 
age velocity V{x,t). The most well-known model, some- 
times called the LWR model, was proposed by Lighthill, 
Whitham, and Richards [L^TTB"]. It is based on the conti- 
nuity equation 

dp{x,t) .s ggQM) / dV(x,t) 

for the density and a speed-density relationship 

V(x,t) = V e {p(x,t)) (2) 

or, alternatively, a "fundamental diagram" Q(x,t) = 
Q e (p(x,t)) for the vehicle flow Q(x,t) — p(x, t)V (x, t). 
Obviously, the LWR model is based on a (hyperbolic) par- 
tial differential equation of first order. A detailed analysis 
is given in Refs. [12JT4]. It is well-known, that it describes 
the generation of shock waves characterized by discontin- 
uous density changes. 

Therefore, in his famous "Requiem for Second-Order 
Fluid Approximations of Traffic Flow" TjT], Carlos Da- 
ganzo correctly notes on page 285 that, "Besides a coarse 
representation of shocks, other deficiencies of the LWR 
theory include its failure to describe platoon diffusion 
properly ... and its inability to explain the instability of 
heavy traffic, which exhibits oscillatory phenomena on the 
order of minutes." However, he also criticizes theoreti- 
cal inconsistencies of alternative models, which, at that 
time, were mainly second-order models containing diffu- 
sion, pressure, or viscosity terms. The Payne- Whitham 
model [T4 j ll5 | fT6]. for example, has a dynamic velocity 
equation of the form 



dV(x,t) 



at 

V 



+ V{x,t) 
dp(x,t) 



dV(x,t) 
dx 



p(x, t) dx 



V e (p(x,t))-V(x,t) (3) 



with 



1 dV c {p) 
2t dp 



1 

27 



dV c (p) 



dp 



> 0. 



(4) 



Here, the term containing v is called anticipation term, 
while the last term is known as relaxation term. V c (p) 
denotes the equilibrium velocity and r the relaxation time. 

Some of the second-order models, including the Payne- 
Whitham model p"5l[T6] , can be derived from car- following 
models by certain approximations. This involves gra- 
dient expansions of non-local, forwardly directed (i.e. 
anisotropic) vehicle interactions [10]. Such approximations 
are problematic, since they lead to terms containing spa- 
tial derivatives, which imply undesired backward interac- 
tion effects as well. The related theoretical inconsistencies 
were elaborated by Daganzo. In the following, we will sum- 
marize his critique by quotes from Ref. [llj (page numbers 
in square brackets): 



1. Lack of anisotropy: "A fluid particle responds to 
stimuli from the front and from behind, but a car is 
an anisotropic particle that mostly responds to frontal 
stimuli" [p. 279]. 

2. Insufficient description of jam fronts: "The width 
of a traffic shock only encompasses a few vehicles", 
while second-order models involving viscosity terms 
would typically imply extended jam fronts [p. 279]. 
Daganzo argues that "the smoothness of the shock is 
inherently unreasonable" [p. 282], because "spacings 
and density must change abruptly whenever the road 
behind is empty" [p. 282]. Based on the analysis of con- 
crete examples, Daganzo further finds that "the cars 
at the end of the queue move back and the behavior 
spreads to the remaining vehicles in the queue ... from 
the back to the front!" [p. 283] . Further on, new arrivals 
of vehicles would "compress a queue from behind" [p. 
283]. 

3. Insufficient representation of acceleration pro- 
cesses and driver characteristics: According to the 
"relaxation" mechanism for the velocity distribution 
assumed in the gas-kinetic traffic model by Prigogine 
et al. |17j , the "desired speed distribution is a property 
of the road and not the drivers, as noted by Paveri- 
Fontana (1975)" [p. 280]. However, "Unlike molecules, 
vehicles have personalities (e.g., aggressive and timid) 
that remain unchanged by motion" [p. 279], and mod- 
els should make sure "that interactions do not change 
the 'personality' (agressive/timid) of any car" [p. 280]. 
Therefore, "a slow car should be virtually unaffected 
by its interaction with faster cars passing it (or queue- 
ing behind it) ..." [p. 280]. 

A further criticism concerns the propagation speeds of 
perturbations in the traffic flow, predicted by second- 
order traffic models, which will be addressed after we have 
replied to the above, well-taken points: 

1. The lack of anisotropy is a consequence of gradient ex- 
pansions and can be avoided by non-local macroscopic 
traffic models [10] , such as the gas-kinetic-based traffic 
model (GKT model) PITTS] . 

2. Non-local traffic models can represent sharp shock 
fronts well, as has been demonstrated for the GKT 
model [20] . They are also capable of avoiding negative 
vehicle velocities, if properly specified TJU]. For exam- 
ple, the speed variance 8 appearing in some macro- 
scopic traffic models, in particular in the "pressure 
term" (see below) must vanish, whenever the average 
velocity V vanishes. This can be reached by a rela- 
tionship of the form 8(p, V) — a(p)V 2 with a suitable, 
density-dependent function a(p) > plfT9] , 

3. The personality of drivers can be represented by multi- 
class traffic models [T9lf21][2~2l . Moreover, the unrealis- 
tic acceleration-behavior implied by Prigogine's gas- 
kinetic traffic model |17] has been overcome by the 
gas-kinetic model by Paveri-Fontana T23] and its gener- 
alizations to different driver- vehicle classes 19,21 . In 
these models, it is not the velocity distribution which 
relaxes to a desired velocity distribution (which would 
imply discontinuous velocity jumps at a certain rate). 
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Rather they describe a continuous adaptation of indi- 
vidual vehicle velocities to their desired speeds. 

Let us now turn to the discussion of the "characteristic 
speeds" . Characteristic speeds relate to the eigenvalues of 
hyperbolic partial differential equations. They determine 
the solutions for given initial and boundary conditions, in 
particular which locations influence the solution at other 
locations at a given time [24]l2"5] (see Appendix |Aj . The 
characteristic speeds are also important for the stability 
of numerical solution schemes for partial differential equa- 
tions [26] . 

What implications does this have for macroscopic traf- 
fic models based on systems of hyperbolic partial differ- 
ential equations with source terms? In his "Requiem for 
second-order fluid approximations of traffic flow" [IT] , Da- 
ganzo argues that "high-order models always exhibit one 
characteristic speed greater than the macroscopic fluid ve- 
locity. ... This is highly undesirable because it means that 
the future conditions of a traffic element are, in part, de- 
termined by what is happening ... BEHIND IT! ... it is a 
manifestation of the erroneous cause and effect relation- 
ship between current and future variables that is at the 
heart of all high-order models" [p. 281]. 

Is this violation of causality a result of crude approxi- 
mations underlying second-order macroscopic traffic mod- 
els? Or could the assumption of circular boundary condi- 
tions explain an influence from behind, even in the case 
where vehicle interactions are exclusively directed to the 
front? Or is the faster characteristic speed related to vehi- 
cle interactions at all? Until today, the problem of charac- 
teristic speeds is puzzling, and it has stimulated many sci- 
entists to develop and investigate improved macroscopic 
traffic models [^[MlllllSniinilSlSllSllSSlISS]- Here, we 
restrict our discussion to the most prominent example: 
In their "Resurrection of 'second order' models of traffic 
flow" [57], Aw and Rascle propose a new model with two 
characteristic speeds, one of which is smaller than and the 
other one equal to V, where V denotes the macroscopic 
vehicle speed. Details are discussed in Sec. 14.11 While, 
without any doubt, such an approach is interesting and 
worth pursuing, we will address the question, whether it 
is necessary to overcome the problem pointed out by Da- 
ganzo. This issue must be analyzed very carefully in order 
to exclude misunderstandings and to avoid jumping to a 
conclusion. To provide a complete chain of arguments, the 
main text of this paper is supplemented by several appen- 
dices. 



3 Linear Instability of Macroscopic Traffic 
Models 

Let us start our analysis with the continuity equation (TTJ) 
for the vehicle density p(x, t) and a macroscopic equation 
for the average velocity V(x, t) of the type derived at the 
end of Sec. 4.4.3 of Ref. [TIT]: Assuming repulsive vehicle 
interactions that depend on the vehicle distance and vehi- 
cle speed, but (for simplicity) not on the relative velocity, 



it reads 



at ox 

l dPjjp, V) dp{x,t) _ l dP 2 (p, V) dV(x,t) 
p dp dx p dV dx 
V Q {p, V)-V(x,t) 



(5) 



Herein, Pi and P2 are contributions to the "traffic pres- 
sure", and V Q (p, V) is the "optimal velocity" function. 

Our stability analysis starts with an initial state of uni- 
form vehicle density p e . The related stationary and homo- 
geneous (i.e. time- and location-independent) solution is 
obtained by setting the partial derivatives d/dt and d/dx 
to zero. In this way, Eq. ([5]) yields the implicit equation 



vM = v { Pe ,vM) 



(G) 



for the equilibrium speed V e (p e ). With this, we can define 
the deviations 



Sp(x,t) = p{x,t) - p e 



and 



5V(x,t) = V(x,t)-V e . 

(7) 

Inserting p (x, t) = p c + 5p(x,t) &ndV(x,t) = V c + 5V(x,t) 
into the continuity equation, performing Taylor approx- 
imations, where necessary, and dropping all non- linear 
terms because of the assumption of small deviations 
8p(x,t)/p c <C 1 and 5V(x,t)/V c <C 1, we end up with 
the following linearized equation: 



d5p(x,t) dSp{x,t) 

f Vo(/0 e 



dt 



dx 



dSV(x,t) 
' Pc dx ■ (8) 



Analogously, the linerarized dynamical equation for the 
average velocity becomes 

dSV(x,t) t 86V(x,t) 



dt dx 
dPi(p e ,V e )dSp(x,t) , dP 2 (p c ,V c )dSV(x,t) 



dp 

dv ( Pe ,v e ; 



dp 

dV (p e ,Ve) 

8V 



dx 
Sp(x,t) 

5V(x,t) - 5V(x,t) 



dV 



dx 



(9) 



The terms on the right-hand side in the first square 
bracket may be considered to describe dispersion and 
interaction effects contributing to the "traffic pressure" , 
while the terms in the second square bracket result from 
the so-called relaxation term, i.e. the adaptation of the av- 
erage velocity V(x, t) to some "optimal velocity" V (p 7 V) 
with a relaxation time r. 

As is shown in Appendix |B] a linear stability analysis 
of Eqs. ([5]) and (JHJ) leads to the characteristic polynomial 



(A) 2 + A 



\np c 



m dP 2 

in dPi 
pe dp 



T \ dV 

I dVo 
t dp 



= 0. 



(10) 
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It has the two solutions (eigenvalues) 
1 in dP 2 



2f 2p e dV 



±y/K±i\$S\ (11) 



with 



-J_ = Ifl-^)>0 
f( Pe ,K) r \ L 8V ] ■ 



(12) 



1 ,<9Pi k 2 fdP 2 s - 



^M = TH-^--^r-£\ , (13) 



±|9f(p e ,K) 



4f 2 dp 4p c 2 

K/3 e rfV'o ft dP 2 



t dp 2p e i dV 
Here, we have used the abbreviations 



(14) 



A = A — kl> and Co = lo — KV c (p c ) . (15) 

As the square root contains a complex number, it is 
difficult to see the sign of the real value A of A. However, 
we may apply the formula 



2 + 3 2 + 5R 



±iy-(ys 2 + 3 2 -5R) , (i6) 

which is derived in Appendix [Cj From this and Eq. , 
we get the following relationship for the real part of the 
eigenvalues A±(p c , k): 



A±(po, k) = Re(A± ( Pc , K )) = -— ±J- + & + 3? 

( 1? ) 

The expression for the imaginary part gives 
- Co±(p c ,K) = Im(A±(p e , «)) 



•^ ± V5^* ! + - VJ - ,, I- 



3.1 Derivation of the Instability Condition 

A transition from stable to unstable behavior, i.e. the 
change from negative to positive values of \±{p e ,K) oc- 
curs only for the eigenvalue \ + (p Q ,n), namely under the 
condition 



\ + { Pe , K ) = -L + ^l (yu 2 +^ + u) =o. (19) 

This implies 
and, therefore, 



1 * °t 2 

+ (21) 



Inserting the above definitions of 9 and 5ft, we eventually 
find 



K 

if 2 
1 

4 



dPi i / ap 2 

9p 4p c 2 \dV 

Kfc&Vo k dP 2 
t dp 2p c r dV 



(22) 



From this and definition (fT2|) , we can derive the following 
condition for the instability threshold: 




PcdVo 
t dp 



1 dP 2 
2 Pc t dV 



(23) 



Assuming the relationships dV (p)/dp < 0, dV /8V < 0, 
and dP 2 /dV < 0, the condition for Re(A + ) > becomes 



dV 



dp 



> 



ldP 1 


H 1 


1 


dP 2 


/ dp ' 


4p c 2 


'We 


dV 



dV 



dV 



(24) 



16f 4 4f 2 4 



We notice that this instability condition is not fulfilled, if 
the average velocity V (p, V) changes little with the den- 
sity p, which is typically the case for small densities and, in 
many models, also for large ones. However, A+(p c , n) may 
be greater than zero at medium densities, where \dV c /dp\ 
is large according to empirical observations. The related 
instability mechanism is based on a reduction of the aver- 
age velocity with increasing density. Due to the continuity 
equation, this tends to cause a further compression (but 
the "traffic pressure" terms P\ and P 2 partially counteract 
this re-inforcement mechanism). 

As a consequence of the inequality (|24|) . we can state 
that the speed-dependence of the traffic pressure term P 2 
and the optimal velocity V Q tends to make traffic flow 
more stable with respect to perturbations. The speed- 
dependence also resolves problems related to the fact 
that dP\jdp may become negative in a certain density 
range. This would imply a negative discriminant of the 
square root, if the negative contribution dP\jdp < was 
not compensated for by (dP 2 / 'dV) 2 j '(4p e 2 ) \1Q\- The case 
dP\jdp < could also cause negative accelerations and 
speeds, particularly at the end of congestion areas, which 
would not be realistic [11] . Again, the second pressure con- 
tribution P 2 can resolve the problem, if properly chosen. 



3.2 Characteristic Speeds, Phase, and Group Velocities 

When neglecting the relaxation term (i.e. in the limit 
t — > oo) , the so-called characteristics may be imagined 
as (parametrized) space-time lines, along which the solu- 
tion of a macroscopic traffic model based on partial dif- 
ferential equations does not change in time. In Appendix 
lAl we derive the characteristics of the linearlized equa- 
tions ((8]) and ([9]). In the following, we will compare the 
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characteristic speeds Cj{p c ) = V c (p c ) + Cj(p c ) given by 
Eq. (|66|) with the phase velocities V e (p c ) + cj±(p Cl k)/ k 
and the group velocities V c (p c ) + duj±{p c , k)/8k result- 
ing from the above linear instability analysis. While the 
phase velocity describes the propagation of a single wave 
mode, the group velocity describes the propagation of a 
wave packet composed of waves with different wave num- 
bers k (see Appendix |D] for details). The group velocity 
is usually considered to represent the speed of informa- 
tion propagation^ Due to dispersion effects, we may have 
du±(p e ,K)/dn ^ &±(p e ,K)/K. 

Let us first study the situation in the limit r —> oo of 
arbitrarily slow adaptation to changed traffic conditions. 
Considering the definitions (fT2|) to (fl"4f . we find 1/t(k) = 
0, |9f(p e ,re)| = 0, and 



sol \ i dPl ^ , ■ - 



dp Ap c 



dV 



(25) 



For 5i < 0, we have V^ 2 + 3 2 
Eqs. (d?!) and $TE§, we obtain 



and, due to 



A± = and 



K 

2pZ 



dP 2 



dV 



in the limit r — ► oo. This implies 

d&±(pe,K) UI±(p e ,K) 



(26) 



Ok 



1 

'2pZ 



dP 2 



dV 



Therefore, group and phase velocity in the limit r — > oo 
are the same. A comparison with Eq. (|6"6"|) shows that they 
also agree with the characteristic speeds. This is expected, 
because of X± =0, which means that the wave amplitudes 
do not grow or decay — they just propagate along the char- 
acteristics. 

For finite values of r, which are typical for real traffic 
flows, the phase and group velocities may be different, 
and they also do not need to agree with the characteristic 
speeds, as we will see below: The group velocities, i.e. the 
propagation speeds of small perturbations, are given by 



du)i(p e ,K) dui(p c ,n) 

C t (p c , K) = = Vc(Pc) + 



dn 

= V c (p c ) + Ci(p e , k) . 



8k 



(28) 



as derived in Appendix [Dl Obviously, there are two group 
velocities C± = V c + c±, which can be determined by 
differentiation of the expression for u>±{p c , k) given in Eq. 



C±(Pc, K) 



1 3P 2 d /I 

T ■ 



2p e dV 1 Ok V 2 



9? 2 + 3 2 - . (29) 



1 A typical example is the modulation of electromagnetic 
waves used to transfer information via radio. 



Considering dP2/dV < and 

1 

2 



5R 2 



= I A± + i ] - )1 



(30) 



which is implied by Eqs. (| 1 T[) and (|18p . we may also write 
dP 2 



C±(Pe 



1 

2p~c 



dV 



d_ 

dn 



A+ + 



-K. (31) 



Taking into account Eq. (|T3|) . this is generally not the 
same as uj±(p c , k)/k, i.e. the phase velocities differ. Inter- 
estingly enough, however, at the stability threshold given 
by A + = 0, we find 



1 



C + iPc > K) = ~2p- c 
1 



dPo 



dV 
dP 2 



dV 




(32) 



At the stability threshold we furthermore have A_ 
— 1/f. Inserting this into Eq. (f3Tj) reveals 



c_ (p c ,K) 



1 

2p~c 
1 

'2p~c 



dP 2 



dV 
dP 2 



d_ 

Ok 



1 

4f2 



5R 



dV 



dp 4p 2 



dP 2 
dV 



(33) 



The same expressions are found for the phase velocities. A 
comparison with Eq. (pB]) shows that they also agree with 
the characteristic speeds. Note that c + is smaller than 
zero. However, we have c_ < (corresponding to charac- 
teristic speeds slower than the average vehicle velocity or 
equal to it) only if 



(34) 



or 




(35) 



4 Discussion 

For the discussion of our results regarding the character- 
istic speeds, let us study two particular models first, the 
Payne model [15lfT6] and the Aw-Rascle model [27] , 



4.1 Characteristic Speeds in the Aw-Rascle Model 

The model proposed by Aw and Rascle [57] corresponds 
to Eqs. H]) and ([5]) with r — > oo, 



dPijp, V) 
dp 



and 



dP 2 ( P , V) 
dV 



- JP (x, ty +1 <o, 

(36) 
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see Ref. [10]. 7 is a positive constant. This implies 1/r = 
0, K(k) = -K 2 {dP 2 /dV) 2 /{4p c 2 ) < and |9f(/s)| = 0. 
Therefore, Eq. (j2"9"]) implies 



C±(p e ,«) = - — 



1 

2p e 
1 

2^ 



ap 2 



dV 
dP 2 



dV 



1 

2pZ 



dP, 



dV 



(37) 



This leads to c + = — 7/5(2;, i) 7 and c_ = 0, correspond- 
ing to the characteristic speeds V — 7/0(2;, £) 7 and V", in 
agreement with Aw's and Rascle's calculations 27 . That 
is, their model does not have a characteristic speed faster 
than the average vehicle speed, which elegantly avoids the 
problem raised by Daganzo [IT] . 

However, is it really necessary to exclude the existence 
of a characteristic speed faster than the vehicle speeds? In 
order to address this problem, we will now study Payne's 
macroscopic traffic model, which has received most of the 
criticism. We do this primarily for the sake of illustra- 
tion, while we are well aware of the weaknesses of this 
model (like the possibility of backward moving vehicles at 
upstream jam fronts for certain initial conditions). There- 
fore, the authors of this paper generally prefer the use of 
non-local macroscopic traffic models [10], but this is not 
the issue to be discussed, here. 



4.2 Payne's Traffic Model 

Payne's macroscopic traffic model [15 1 116 ] has a solely 
density-dependent optimal velocity 



V ( P , V) = V e (p) 



and the pressure gradients 



dP^V) _ 1 
dp 2r 



dV c {p) 



dp 



> 0. 



dP 2 (p, V) 



dV 



(38) 



= 0. (39) 



This simplifies the instability condition 
and we get 



Pc 



dV c (p c 



dp 



> 



1 



2p c r 



considerably, 



(40) 



Traffic flow becomes unstable, if the equilibrium velocity 
V (p) decreases too rapidly with an increase in the density 
p, and greater relaxation times r tend to imply larger 
instability regimes. For the characteristic speeds at the 
instability threshold, with p c \dV c /dp\ = l/(2p c r) we find 



c±0b) = T\ -5— 




TPc 



dV c (p c ) 



dp 



(41) 

Clearly, c_(p) is non-negative, i.e. the related characteris- 
tic speed V c (p) + c-(p) tends to be larger than the av- 
erage vehicle speed V e (p). Nevertheless, by demanding 
V c (p) + c_(p) < V°, e.g. by assuming a linear speed- 
density function 



(42) 



V c (p) = V" 1 



Pja 



one could still reach that the characteristic speed V c (p) + 
c_ (p) lies within the variability of the vehicle speeds. In 
fact, we have c± — whenever the vehicle speed cannot 
vary, namely at density zero and at maximum density, 
where p c \dV c (p c ) / dp\ = 0. However, do we need to im- 
pose such conditions on the characteristic speed and the 
speed-density relationship? This shall be addressed in the 
following and in Sec. [5] 

In connection with this question, it is interesting to 
note that, according to Eqs. (f33j) and (|41[) . the group ve- 
locity c + corresponding to the solution with the unstable 
eigenvalue A + is negative with respect to the average ve- 
locity V c . In contrast, propagation at the positive speed 
c_ with respect to the average velocity V e is related with 
an eigenmode that decays quickly, basically at the rate at 
which the vehicle speeds adjust. Therefore, the forwardly 
propagating mode cannot emerge by itself. It could only be 
produced by a particular specification of the initial condi- 
tion, enforcing a finite amplitude of the forwardly moving 
mode. We will come back to this in Sec. [5] 

It is noteworthy that already Whitham performed a 
thorough analysis of the speeds characterizing the traffic 
dynamics in what is known as the Payne model today (see 
Ref. rj_2], Chaps. 3 and 10). He showed that the linearized 
partial differential equations |(8]) and ([9]), when specified 
in accordance with Eqs. I|38p and ([3"9"]> . can be cast into 
the equation 

dSp(x,t) ,( V (,, dV c (p)\ dSp(x,t) 



dt 



dp 



Ox 



- r (^ +[Fe(P)+C+(P)] ^ 
x(§- t + [V c (p) + c-(p)]^)5p(x,t). 



(43) 



Whitham was perfectly aware of the fact that the char- 
acteristic speed V c (p) + c_(p) was faster than the average 
vehicle velocity V c (p), but not at all worried about this. 
His perception was that all three velocities were meaning- 
ful, and that the kinematic speed V c (p) + pdV c /dp would 
dominate in the limit of small values of r (which implies 
stable vehicle flows). However, the open problem is still, 
how a characteristic speed V c (p) + c-(p) > V c (p) can be 
interpreted, without violating causality. 



4.3 Characteristic Speeds vs. Vehicle Speeds 

In physical systems, it is not necessarily surprising to find 
characteristic speeds faster than the average speed. Let 
us illustrate this for the example of sound propagation. 
In one spatial dimension, this is described by the continu- 
ity equation ([T]) in combination with the one-dimensional 
velocity equation 



dV(x,t) 
di 



V(x,t) 



dV(x,t) 
dx 



ldV(p) 
P 



dx 



(44) 



These so-called Euler equations [37] can be considered 
to model frictionless fluid or gas flows in one dimension. 
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Compared to the velocity equation ([5|), we have dropped 
the relaxation term [V (p) — V]/t. Therefore, we do not 
have an equilibrium velocity-density relation V e (p), now. 

In order to determine the solution of the above equa- 
tions, one can derive linearized equations for the case of 
small deviations Sp(x,t) — p(x,t) — p c and SV(x,t) — 
V(x, t) — V e from the stationary and homogeneous solu- 
tion p(x,t) = p c and V(x, t) = V e = 0. The quantity p c 
corresponds to the average density of the fluid or gas. 

Inserting (J7JI into Eqs. |T|) and (14"4"1) and neglecting non- 
linear terms in the small deviations Sp, SV results in 

dSp(x,t) dSp(x,t) dSV(x,t) 

V e - = -p e — (45) 



dt 



Ox 



dx 



and 



dSV(x,t) , „ dSV(x,t) 1 dV( Pc )dSp(x,t) 

Vo = : ■ (46) 



Of 



dx 



p c dp dx 



Konhauser [39], and Lee et al. [40] . Their pressure func- 
tions are also given by the formula "density times veloc- 
ity variance" . Therefore, the faster characteristic speed of 
these macroscopic traffic models is expected to lie within 
the range of individual vehicle speeds o 

As we have seen above, the situation is generally differ- 
ent for Payne's model. However, it is illustrative to note 
that V a (p) + c + (p) may become negative, even when all 
vehicles move forward. That is, it is possible to have char- 
acteristic speeds outside of the range of vehicle speeds: 
According to Eqs. (|4"Tj) and (fT5|) . the slower characteristic 
speed at the instability threshold is 



V e (p) + cjp) = Vc( P ) - p 



dVc{p) 



dp 



Ve(p)+P 



dV e ( P ) dQc(p) 



dp 



dp 



(52) 



Considering V c — 0, deriving Eq. (|45j) with respect to t, 
and Eq. (|46[) with respect to x yields 



d 2 Sp(x,t) d 2 SV(x,t) n 

+ Pc =0 



dt 2 



dt dx 



and 



d 2 8V{x,t) _ 1 dV{ Po ) d 2 Sp{x, t) 



(47) 



(48) 



dx dt p c dp dx 2 

Inserting Eq. ijlgjl into Eq. ([4"T]) finally gives the so-called 
wave equation 



d 2 Sp(x,t) g2 d 2 5p(x,t) = Q 



dt 2 



dx 2 



(49) 



which is well-known from one-dimensional sound propa- 
gation. The constant 



UVipc) 



dp 



(50) 



corresponds to the speed of sound. In order to determine 
the spatio-temporal solution of Eq. (|4"9"|) , we rewrite this 
equation, inspired by the relationship (a 2 — b 2 ) = (a + 
b)(a-b): 



(-- 


- 9 \ 


(~ 




\dt 


^Tx) 


\dt 


-*Tx) 



Sp(x,t) = 0. 



(51) 



According to this equation, perturbations propagate back- 
ward and forward at the speed ±c, although the average 
speed is V = 0. However, for gases we may assume an ap- 
proximate pressure law of the form V = p6o [37| , where 9q 
is the velocity variance of gas molecules. Hence, the speed 
of sound is given by c = ^/0q, i.e. by the standard devia- 
tion of velocities. As a consequence, the speed of sound can 
actually be propagated by the mobility of gas molecules. 

In a similar way, we can understand characteristic 
speeds faster than the average vehicle speed in the macro- 
scopic model of Phillips [38] or Kiihne [8], Kerner and 



Since Q c (p) = pVc(p) represents the "fundamental dia- 
gram", dQ e (p)/dp describes the negative speed of kine- 
matic waves in the congested regime [13]. This does not 
constitute any theoretical inconsistency, even if Ve{p e ) + 
c + (p) < 0. In fact, we all know situations involving neg- 
ative group velocities from dissolving congestion fronts, 
e.g. when a traffic light turns green: There, the negative 
propagation speed just results from the fact that the con- 
gestion front moves backward, whenever vehicles leave a 
congested area with some delay. Hence, the negative char- 
acteristic speed does not describe the speed of cars. It 
reflects the propagation of gaps rather than vehicles. 

Therefore, could we have a similar mechanism that 
generates characteristic speeds faster than the vehicle 
speeds? If vehicles would react to their leaders with a neg- 
ative delay, this would in fact be the case, but it would 
violate causality. Therefore, all possible explanations for 
characteristic speeds faster than the vehicle speeds con- 
sidered so far have failed to resolve the problem. However, 
the problem may still be a result of the approximations 
underlying second-order macroscopic traffic models. As we 
have indicated before, the gradient expansion required to 
derive them implies some degree of backward interactions. 
Therefore, it is conceivable that following vehicles would 
cause their leaders to accelerate, even beyond their desired 
speed V°. 

If this would be the explanation of a characteristic 
speed faster than the average speed V or free speed V°, 
we should not observe it in microscopic traffic models 
with forward interactions only. Therefore, we will now de- 
termine the characteristic speeds of the optimal velocity 
model [9 . This car-following is chosen, because the Payne 
model can be considered as a macroscopic approximation 
of it (see [TU] and references therein). Besides, we will com- 
pare the instability conditions of both models. 



2 Note that the existence of perturbations in the traffic flow 
always implies a variation of the vehicle speeds. 
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5 Linear Instability and Characteristic Speeds 
of the Optimal Velocity Model 

We have seen that macroscopic traffic models behave un- 
stable with respect to small perturbations in a certain den- 
sity range, where the average velocity changes too rapidly 
with the density. The same is true for many car-following 
models. As an example, we will shortly discuss the dy- 
namic behavior of the optimal velocity model. While its 
stability has been already studied in the past [9], we will 
focus here on the characteristic speeds, in order to show 
that characteristic speeds greater than the average veloc- 
ity are not an artifact of macroscopic traffic models. 

According to the optimal velocity model, the change 
of the speed Vi (t) of vehicle i is given by 



dvj _ v (dj(t)) - Vj(t) 
~dt ~ t 



(53) 



and the temporal change of the distance di(i) = Xi-i(t) 
Xi (i) to the leading vehicle i — 1 is determined by 



ddl M 



Vi{t) 



(54) 



In the above equations, the distance-dependent function 
v Q (di) is called the optimal velocity function and r is again 
the relaxation time for adjustments of the speed. 

Appendix IE1 sketches the linear stability analysis of the 
optimal velocity model. In the following, we will focus on 
the analysis of the group velocity c± with respect to the 
average velocity v Q (d c ), i.e. the velocity at which perturba- 
tions are expected to propagate. Relative to the average 
motion of vehicles with speed v c (d c ), the characteristic 
speeds are 



c±(d e , k) 



dCb±(d c , k) 
dn 

-LJL 
2^dk V 2 



L dCj±(d e ,k) 
2^ dk 



(55) 



This can be derived analogously to Eq. ([29)1 . using Eq. 
fT6]) and n = 2nk/L. According to Eq. (f3Tj) and due to 
the series expansion cos(a;) w 1 — x 2 /2, at the instabil- 
ity threshold with A + = and dv (d e ) / dd = l/(2r), we 
obtain with Eq. ([TOST) 



T ^JlJ^ [1 _ co8{2nk/N)] 



2tt dk 



dd 



L d jl dv a {d c ) 1 /27rfeV 
27r ~dk y r dd 2 [~N~ ) 



N V 2r dd 



Td e 



{ dv (d c ) 
V dd 



Td ( 



2r dd 

dv (d c ) 
dd ■ 



(57) 



It is remarkable that the group velocity of the optimal ve- 
locity model can again exceed the average vehicle velocity 
v (d e ), namely by an amount c_(d e ) = d c dv (d e )/d e > 0. 
Moreover, it can be shown that the instability thresholds 
and the related characteristic speeds are the same as for 
the Payne model (see Appendix IF]) . This confirms that the 
Payne model may be viewed as macroscopic approxima- 
tion of the optimal velocity model (see [TU] and references 
therein). In view of these results, it is hard to argue that a 
characteristic speed faster than the vehicle speeds consti- 
tutes primarily a theoretical inconsistency of certain kinds 
of macroscopic traffic models. Quite unexpectedly, it also 
occurs for microscopic traffic models that, according to 
computer simulations, behave reasonably well. 

Therefore, the approximations underlying the Payne 
model cannot be the problem for the existence of a char- 
acteristic speed faster than the vehicle speeds. However, 
it is interesting to note that the larger group velocity 
v (d e ) + C-(d e ) is related to a negative real part A of 
the eigenvalue A_. According to Eq. (|29]) , the fast char- 
acteristic speed V e (p c ) + c_(p c ) of macroscopic second- 
order models is related to a negative eigenvalue A_(p e ) as 
well, see Eq. (|17j) . Therefore, the related eigenmode decays 
quickly, and it will be hard to observe in reality. In par- 
ticular, the faster propagating mode may not emerge by 
itself. A closer analysis shows that both, for the optimal 
velocity model and the Payne model, A_ is of the order 
— 1/r, i.e. related to the relaxation time r of vehicles. We 
will see that this observation is highly relevant for under- 
standing perturbations that move faster than the vehicles 
do. 

After all, does the fast characteristic speed really con- 
stitute a theoretical inconsistency? Not so, if we can find 
initial conditions, for which a following car accelerates or 
decelerates earlier than the leading car does, although the 
leader does not react to the follower. In fact, such initial 
condition can be constructed: Figure Q] shows the result of 
a computer simulation with vehicles on a circular road 
of length L. We assume that all vehicles have the distance 
d = d c = L/N initially. Moreover, all vehicles, with the 
exception of 10 subsequent vehicles, are assumed to have 
the initial speed v (d e ). Furthermore, the speed of the last 
of the 10 vehicles is set to (or v°), the speed of the first 
one to v (d c ). The speeds of the vehicles in between are 
determined by linear interpolation. For this scenario, it is 
quite natural that the last of the 10 vehicles accelerates 
(or decelerates) first, since it experiences the largest devi- 
ation of its actual velocity Vi(0) from the optimal velocity 
v (d e ). However, as this earlier acceleration (or decelera- 
tion) is not interaction-induced, it does not violate causal- 
ity. The large characteristic speed in macroscopic traffic 
models can be understood in a similar way. 



6 Summary, Conclusions, and Outlook 

In this paper, we have started with a discussion of Da- 
ganzo's sharp criticism of second-order macroscopic traffic 
flow models [TT] . We have argued that most of the deficien- 
cies identified by Daganzo were fully justified, but could 
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Fig. 1. Simulation result of the optimal velocity model with 
v a (d) = u°{tanh[(d - l)/s - 1.2] + tanh(f .2)}/2, v° = 
f f5 km/h, so = 50 m, and I = 4 m. We have chosen a particu- 
lar initial condition, where all vehicles started with a distance 
d c = 200 m to their respective leader, but some vehicles i had a 
speed Vi(0) < v a (d c ) in the beginning. As a consequence, these 
vehicles adjusted their speeds to the optimal velocity. The rel- 
evant point here is that followers reach the optimal velocity 
(or certain fractions of it) earlier than their respective leaders. 
That is, for the particular initial condition chosen here, the 
perturbation in the speeds propagates faster than the vehicle 
speeds. This effect, however, does not violate causality, as the 
earlier acceleration of upstream cars is not triggered by inter- 
actions with followers — it just results from the relaxation term. 
Therefore, the perturbation disappears on a time scale that is 
determined by the relaxation time r = f s, as predicted by the 
real part of the eigenvalue A_, see Eq. (|f04[) . The relaxation 
takes longer for larger values of r. In the limit r — > oo, the 
perturbation does not decay anymore, but according to Eq. 
(|f 04[) . we then have c± — > 0. Therefore, despite its fast speed, 
the perturbation did not overtake the first car upstream of the 
initial perturbation in our simulations, when the parameters 
were chosen in a way that avoided accidents. This confirms the 
validity of the causality principle. 



be overcome in the course of time by improved macro- 
scopic traffic models, particularly by non-local multi-class 
models. However, the issue of characteristic speeds faster 
than the average vehicle speed was still an open, contro- 
versial problem, as it seems to violate causality. In order 
to study it, we have performed a linear instability analy- 
sis of a generalized macroscopic traffic model, which took 
into account speed-dependencies of the optimal velocity 
and the traffic pressure terms. Such speed-dependencies 
occur, for example, in Aw's and Rascle's model [57]. They 
result when realistic vehicle interactions are considered, 
and when the possibility of accidents and negative vehicle 
speeds shall be avoided [4llfTU] . Requirements for reason- 
able models seem to be 



dVojp, V) 

dp 



<o, 



dVo(p,V) 
dV 



< 0, 



dP2 (P,V) 
dV 



<0, 
(58) 



dp 



1 



dP 2 (p, V) 
dV 



> 0. 



(59) 



These conditions are, for example, fulfilled by the gas- 
kinetic-based traffic model (GKT model), see Ref. [4"3"] . 

Our main attention was dedicated to the characteris- 
tic speeds (or group velocities) rather than the instability 
thresholds. In the following, we summarize the main re- 
sults: 

1. While the characteristic speeds may generally differ 
from the group and the phase velocities, in the limit 
t — ► oo of a vanishing source (relaxation) term, they 
are all the same. Therefore, using a different definition 
of propagation speeds does not resolve the problem of 
characteristic speeds faster than the (average or max- 
imum) vehicle speed. 

2. Velocity-dependent pressure terms tend to reduce the 
characteristic speeds, see Eq. (|3Tj) . This is best illus- 
trated by Aw's and Rascle's model, where the fast 
characteristic agrees with the average vehicle speed. 

3. Most macroscopic traffic models have a characteristic 
speed faster than the average velocity, but it may still 
be within the variability of the vehicle speeds, see Eq. 
PI) and Sec. E2 

4. In some models like the Payne model, the characteris- 
tic speeds can move slower than the slowest vehicle and 
faster than the fastest vehicle. The first case is related 
to delayed acceleration maneuvers at jam fronts and 
related to gap propagation during jam dissolution, but 
the second case remained a mystery for a long time. 

5. The faster characteristic speed is related with a negative 
real part of the eigenvalue. This causes a quick decay 
of the corresponding eigenmode, basically at the rate, 
at which the vehicle speed is adjusted. Therefore, this 
eigenmode will not emerge by itself (see Sec. I3.2[) . 

6. If the faster characteristic speed were a result of inter- 
actions with following vehicles in a circular road ge- 
ometry (where following vehicles influence the down- 
stream flow as well), the fast eigenmode should decay 
with the length L of the circular road, not with the 
relaxation time r. Therefore, periodic boundary con- 
ditions cannot be responsible for a characteristic speed 
faster than the vehicle speeds. This has also been ver- 
ified with simulations^ 

7. A characteristic speed faster than the vehicle speeds 
cannot be explained as a result of the approxima- 
tions underlying macroscopic second-order models, as 
it is also found for microscopic car-following models, 
in which vehicle interactions are forwardly directed 
and velocities are restricted to a range between zero 
and some maximum speed. For the macroscopic Payne 
model and the optimal velocity model, we have shown 
a correspondence not only of the instability thresholds, 
but also of formulas for the group velocities (see Ap- 
pendix |FJ . 



A Simulations for open boundary conditions basically yield 
the same results as for periodic boundary conditions, given 
the system (in terms of the road length L) is sufficiently large. 
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8. Assuming particular initial conditions, characteristic 
speeds faster than the average vehicle speed could be 
demonstrated to exist in computer simulations, where 
followers accelerate (or decelerate) before their leaders 
do (see Fig. [TJ. As these acceleration (or deceleration) 
processes are induced by artificial initial perturbations 
rather than by vehicle interactions, this does not imply 
a violation of causality. 

Given these findings, we conclude that characteristic 
speeds faster than the average speed of vehicles do not con- 
stitute a theoretical inconsistency of traffic models and do 
not need to be "healed" by particularly constructed traffic 
models|3 From our point of view, the problem is that char- 
acteristic speeds are hard to imagine. In fact, there is no 
direct correspondence to particle or vehicle velocities (see 
Sec. [O] and Appendix [D]). The group velocity is nothing 
more than a matter of phase relations between oscillations 
of successive vehicles in an eigenmode, and the interpre- 
tation as speed of information transmission is sometimes 
misleading. 
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and prepared the figure. 
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A Hyperbolic Sets of Partial Differential 
Equations and Characteristic Speeds 

Let us rewrite Eqs. ([5]) and © in the form of a system of 
linear partial differential equations. With 



S(6p,SV) = - 



dV ( Pe ,V e ) 
dp 

dV ( Pc ,V c ] 



8p{x,t) 



dV 



SV(x,t)-SV(x,t) 



(60) 




An A i2 \ q I 5p(x,t) ' 
An A 22 ) 9x \SV{x,t) 




V c { Pc ) 



1 dP 2 (Po,Va) 

dV 

(62) 



4 Of course, this does not speak against models of the Aw- 
Rascle type. 



As will be shown below, the solution of this system of par- 
tial differential equations is given by the initial condition 
Sp(x, 0) and SV(x, 0). The solution procedure consists ba- 
sically of two steps: On the one hand, we must determine 
the so-called characteristics, and on the other hand, we 
must solve a set of ordinary differential equations to find 
the solutions along them (see Ref. |42] and footnote 3): 
With u(x, t) = (5p(x, t), SV(x, t))' and S = (0, S)' (where 
the prime indicates a transposed, i.e. a column vector), we 
can rewrite Eq. ([6Tj) as 



Bu(x,t). (63) 



du(x, t) | ^ du(x, t) 



dt dx 
The source term can be rewritten as S = B_u(x,t) with 



B = 



Bn B12 
B 2 i B22 , 





1 dV (p B ,V B ) 
r dp 







r ( dv ( Po y c ) _ A 
r \ dv 1 J 



(64) 

Now, let Cj denote the eigenvalues of the matrix A. The 
values of Cj — V c (p c ) + Cj satisfying det(.A — Cjl) = are 
given by the characteristic polynomial 



2_c L m ^9P L _ a 

j Pc dv d P < 



(65) 



which results in 



1 dP 2 
2p c dV 



± 



1 



Ap c 2 V dV 



dp 



Furthermore, let Zj be the eigenvectors related with the 
eigenvalues Cj = V c + Cj, i.e. 



Azj — CjZj 



(67) 



Finally, let R = (Rij ) be the matrix containing the eigen- 
vectors Zj as their jth column, and y(x,t) = R~ 1 u(x 1 t) 



or u(x,t) = Ry(x,t). Then, inserting this into Eq. (J63J) 
and multiplying the result with the inverse matrix RT of 
R yields 



9 Vj( x ,t) , c dyj(x,t) 



dt 



dx 



= = {Rr l BRy) 3 . (68) 



For S = (corresponding to the limiting case r — > 00), 
we have 

y j {x ) t)=y j (x-C j t,0), (69) 

which means that the solution does not change in time 
along the characteristics Xj (t) = Cjt. The quantities Cj 
are called the characteristic speedsO If u(x, 0) is the ini- 
tial condition, the solution of the set of partial differential 



5 The idea behind the characteristics is to introduce a param- 
eterization t(si,S2), x(si,S2), which is defined by dt/dsj = 1 
and dx/dsj — Cj. Then, one can rewrite Eq. (|68[) as 



dyj _ dyj(x,t) dt dyj(x,t) dx 



dsi 



dt ds 



+ 



dx 



di- = cr^flv). 



In the generalized coordinates si and S2, the partial differential 
equations in x and t we were starting with, turn into ordinary 
differential equations. These are much easier to solve. 



Dirk Helbing and Anders Johansson: On the Controversy around Second-Order Traffic Flow Models 



11 



equations is 



m(x, t) 



Cjt, 0) 



(70) 



with y(x,0) = RT 1 u(x, 0)0 Therefore, the spatio- 
temporal solution u(x, t) is fully determined by the initial 
condition. In other words, the future state of the system is 
given by its previous state, and the principle of causality 
should be valid. 

B Stability Analysis for Macroscopic Traffic 
Models 

In order to understand the dynamics of traffic flows, it is 
important to find out whether and under what conditions 
variations in the traffic flow can grow and eventually cause 
traffic congestion. For this, it is useful to make the solution 
ansatz 

Sp(x, t) = S Po exp (ikx + (A - iu)t) = 5 Po e xt e ^ x -^ , 

SV(x, t) = SV exp (ikx + (A — ico)t) = SV e xt e ^ KX '^ . 

(71) 

Because of exp(inx) = cos(kx) + isin(/ix) (see Appendix 
IC[) . ansatz ([7T|) assumes that the perturbation of the sta- 
tionary and homogeneous traffic situation can be repre- 
sented as a periodic function with the wave number k and 
wavelength 2tt/k. The wave frequency of Eq. (|7ip is to, 
while Spo exp(Ai) and SVq exp(Ai) are the amplitudes at 
time t. That is, if the "growth rate" A is greater than zero, 
even small perturbations will eventually grow, which can 
give rise to "phantom traffic jams". For A < 0, however, 
the initial perturbation will be damped out and the sta- 
tionary and homogeneous solutions will be re-established, 
i.e. it is stable with respect to small perturbations. 

Below we will see that, for each specification of k and 
the average density p e , there exist two solutions I £ {+, — } 
with the frequencies u>i(k) and the growth rates A/(k). All 
the corresponding specifications of ansatz (fTTj) are solu- 
tions of the linearized partial differential equations. The 
same applies to their superpositions. The general solution 
for an arbitrary initial perturbation is of the form 



Sp(x, t) = / dn Sp (n) exp (ikx + [A;(k) — iw/(K)]i) 

'£{+,-} 

5V(x, t) = Jdn SV l (k) exp (ikx + [A/ (k) - \ui («;)] i) 



(72) 



In order to find the possible K-dependent wave numbers 
u! and growth rates A, we insert ansatz (|71|) into the lin- 
earized macroscopic traffic equations ([8]) and ([9]) and use 



6 Note that formulas (|69[) and (|70|) only apply to the limiting 
case r — » oo, where the relaxation term of the macroscopic 
traffic model vanishes. 



the relationship i 2 = —1. The result can represented as an 
eigenvalue problem: 



'Mn Mi 2 \ /Spo' 
Mai M 22 J \SV , 



! (0 



where 



Mn - -A, 
Mi 2 = -inp c , 

ik dPi 1 dV a 



M 2 i = - 
M 22 = -A 



p c dp t dp 

in 8P 2 1 dV 1 
~p~ c ~dV + t~8V ~ t 



(73) 

(74) 
(75) 

(76) 
(77) 



and 



A = A — iu> with u> = u> — nV c (p c ) . (78) 



Equation (|73|) is fulfilled only for certain values of A(k), 
the so-called "eigenvalues" . These depend on the average 
density p e and solve the characteristic polynomial of sec- 
ond order in A, which is obtained by determining the de- 
terminant 



det(M) = M11M22 - M21M12 



(79) 



of the matrix M and requiring that it becomes zero. The 
corresponding characteristic polynomial is given by Eq. 

(USD- 



C Derivation of Formula (19) 

Remember that a complex number 

z = + i3? = re i<p = r cos(^) + ir $m(ip) (80) 

can be represented in two-dimensional space with coordi- 
nates 3? = Re(z) = rcos(iy9) and 3 = Im(z) = rs'm(ip), 
respectively, called the real part and the imaginary part. 
The absolute value is given as 



r= y/W + 3 2 = y/($t + i90(3t - 13) = sfzJ, = \z\, 

(81) 

where z = 5i — i3 = re ltp is the conjugate complex num- 
ber. The angle ip is determined by 



tan(y) 



sin(y) 
cos((p) 



3 



Im(z) 
Re{z) ' 



(82) 



and the exponential functions is defined as for real num- 
bers by the infinite series expansion 



00 / 
z 



exp(z) = c z = J2j\ 



1=0 



(83) 



where l\ = I ■ (I — 1) . . . 2 • 1. Therefore, the relationships 
for exponential functions apply also to the case of complex 
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numbers, i.e. the product of two complex numbers z\ = 
0?i + iSJi = rie' Vl and z 2 = $t 2 + i^5 2 — T2e 1(p2 is given by 

Zl z 2 = (SR!^ 2 - Si9f 2 ) + i(5Ri3 2 + Si5R 2 ) 



rie I¥,1 r 2 e lv>2 



rir 2 e 



i(ipi+ifi2 



= r x r 2 cos(ipi + tp 2 ) + \r x r 2 sin(</Ji + ip 2 ) ■ (84) 

As the real and imaginary part are linearly independent of 
each other, this implies SJii^ — c ^i c ^2 = T\r 2 cos(</?i + ip 2 ) 
and 3fJi32 + 3i5i2 = rir2sin(</?i + ip 2 ). The inverse of a 
complex number is given by 



1 __ 1 e" 1 ^ 
z re llp r 



The imaginary unit i has the property i 2 
therefore, be written as i = \J — 1 = e 17r ' 2 . 
The square of complex numbers 



(85) 

-f and may, 



re 



r [ cos(</j) ± i sin((/j)] 



can, on the one hand, be written as 

cos 2 (( / s) ± 2icos(< / 9) sin^) — sin 2 (iy9) 



2 2 

z = r 



(86) 



(87) 



On the other hand, using the well-known law e Xl ■ e X2 = 
e x x +x 2 f or f-jjg exponential function, we find the alternative 
representation 

z 2 = r 2 (e ±il ^) 2 = r 2 e ±i2 ^ = r 2 [caa(2tp)±iwa.(2cp)] . (88) 

Comparing the real parts and using the trigonometric re- 
lationship sin 2 (a;) + cos 2 (2;) = 1, we find 

cos(2(^) = l-2sin 2 ((^) = 1-2 [l-cos 2 (<ys)] = 2cos 2 (<^)-l , 

(89) 

from which we can derive the trigonmetric formulas 



sim%/2) = -[1-cosfoO] 



and 



cos 



: (^/2) = -[l + cosM] 



(90) 



(91) 



Therefore, the square root of a complex number is given 

by 

y/z = Vre ±itp/2 = Vr[cos(v/2) ±isin(y>/2)] 



i [r + r cos(v?)] ±iJ-[r-r cos(<p)] . (92) 



Considering = rcos(</?), 3 = rsin(^), and 5R 2 + 3 2 
we end up with the desired equation 



i (V a? 2 + 3? 2 + 3?) ±i y 5 (V a? 2 + 3? 2 - 

(93) 



D Meaning of the Group Velocity 

Let us start with the representation ((72)) of the general 
solution of the linearized system of equations, focussing 
(for simplicity) on the case A/ (ft) = and assuming a 
"Gaussian wave packet" with 



Sp l (n) 



e -( K -K ) 2 /(2<?) 



(94) 



Via the linear Taylor approximation lui(k) — uji(kq) + 
Ci An with Ci = duji(no) I dn and An = (k — Kq), from 
Eq. (72]) we get 



6p(x,t) 

E 
E ■ 
E ■ 



-( K - Ko ) 2 /(2e) 
V2^0 



j[K x— O)((K )i] 



~(zl K ) 2 /(20) 
V2^0 

e -[ziK-ie( a: -c ! t)] 2 /(2e) 



=1 



x e 



-0(x-C,t) 2 /2 



- E 

/£{ + .-} 



e i[ Ko a;-w i (Ko)t] e -e(2;-C ! t) 2 /2 



(95) 



While the single waves of frequency loi(k) move with the 
"phase velocity" x/t = ui(k)/k, it turns out that their 
superposition behaves like a wave with frequency loi(kq) 
and speed x/t = u>i{kq)/kq. However, the wave packet 
or, more exactly speaking, its amplitude e^O-Cit) /2 j g 
moving with the group velocity x/t = Ci = duji(n) / dn. 
Note that the case Ci > lui(ko)/ko, in which the group 
velocity is greater than the phase velocity (wave velocity), 
is possible. It is called "anomalous dispersion" . 



E Linear Stability Analysis of the Optimal 
Velocity Model 

For a linear stability analysis of the optimal velocity 
model, we imagine the situation of N vehicles i distributed 
over a circular road of length L. This allows us to assume 
periodic boundary conditions. The stationary solution for 
this case is given by dvi/dt = and ddi/dt = 0, which 
implies 

di(t) = d c = L/N = const. 
Vi-i(t) = Vi(t) = v (d e ) = const. (96) 

We are now interested how the deviations from this solu- 
tion, i.e. the variables 



Sdi(t) = di(t) - d c , 
6vi(t) = Vi(t) - v Q (d e ) 



(97) 
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develop in time, assuming that the initial deviations are 
small, i.e. <5<fj(0) <C d c and Svi(0) -C v c (d c ). For this, we 
linearize the model equations (|53|) and ([54]) around the 
stationary and homogeneous solution. This results in 



dSvi(t) 1 f dv (d c ) 

dt 
dSdi(t) 

dt 



dd 

5vi-i(t) - 6vi(t) 



8di{t) - 6vi(t) , 



(98) 



For the analysis of stability, we use the solution ansatz 

Svj(t) = 5v ^ k/N+lt = 5v Q e^ L ' N+lt , 

Sdj (t) = Sd e ^/N+~xt = Mq e xj K L/N+\t _ ( 99 ) 

where n = 2nk/L is the so-called wave number, which is 
inversely proportional to the wave length 2it/k = L/k. 
Note that, due to the assumed periodic boundary condi- 
tions, possible wavelength are fractions L/k of the length 
L or the circular road. The shortest wave length is given 
by the average vehicle distance d c = L/N, i.e. k 6 
{1,2,..., N}. Summing up the functions ([9"9"|) over these 
values of k results in the Fourier representation of Svj(t) 
and Sdj(t): 



N 



Svj(t) = Sv^e 1 



2-n-jk/N+Xt 



k=l 
N 



ddj(t) = j2 dd ^' 



2-njk/N+\t 



(100) 



fc=l 



The parameters 5vk and Sdk are determined by the initial 
conditions of all vehicles j. A = A — iui are the so-called 
eigenvalues, whose real part A describes an exponential 
growth (if A > 0) or decay (if A < 0), and whose imag- 
inary part Co reflects oscillation frequencies. Sdg and 5vq 
denote oscillation amplitudes. Inserting this into (|98|) and 
dividing by e l2?T jk/N+\t ^ we g na Uy obtain 



\Sv = - ( —^f~ Sd o -dv ) . 



(101) 



\Sd = 5v e~ i27lk/N - Sv = <^ (c~ i2 " fc/JV - l) . (102) 



Multiplying Eq. (fTUTj) with A and inserting Eq. ([i"02jl for 

X5d(j in the square brackets gives, after division by Svq, 
the characteristic polynomial in the eigenvalues A, namely 



A 2 + I A - 1 dv °( d ° 



dd 



(e" 12 ^-!) =0. 



(103) 



The solutions A(d e , k) of this polynomial are the eigenval- 
ues. They read 



k/N _ I 



(104) 

Again, the square root contains a complex number, which 
makes it difficult to see the sign of the real value A± of A±. 



However, considering e ±IV = cos(c/?)±isin(<£>) and defining 
the real part 



3? = 



1 

4T 2 " 



1 dv (d e 



dd 



[1 - cos(2tt£;/A0] (105) 



of the expression under the root and its imaginary part 



is 



sin(27rfc/iV) dv Q (d e ) 



(106) 



t dd 

we can again apply the useful formula (|16[) . From this we 
can conclude that A = Rc(A) = if 



1 



16t 4 



4t 2 + 4 ' 



(107) 



see Eq. ([21]) • Inserting Eqs. ([1051) and ([TOB]) . we find 



sin 2 (27rfc/Ar) fdv (d) 
4^2 { dd 



1 dv a (d) 
4r 3 dd 



[1- 



which finally results in the condition 



dv {d c ) 
dd 



1 



cos(27rfc/AQ fe=0 J_ 
2r 



rshr(27rfc/A0 



i{2nk/N)] , 
(108) 

(109) 



The limit 2nk/N — > follows from cos^) w 1 — ip 2 /2 and 
sm(ip) w (f in the limit of small wave numbers k — 2nk/L, 
i.e. large wave lengths 27t/k = L/k. 

It can be demonstrated by numerical analyses that 



dv {d c ) 
dd 



> 



2t 



(110) 



constitutes the instability condition of the optimal veloc- 
ity model ([53]) 9J. In other words, if the velocity changes 
too strongly with the distance, small variations of the vehi- 
cle distance or speed will grow and finally cause emergent 
waves, i.e. the formation of one or several traffic jams. 
Since the origin of such a breakdown can be infinitesi- 
mally small, these traffic jams seem to have no origin. 
In such situations, one speaks of "phantom traffic jams" . 
A closer analysis for realistic speed-distance relationships 
v (d) shows that traffic tends to be unstable at medium 
densities p — 1/d, while it tends to be stable at small and 
large densities (where the speed does not change much 
with a variation in the distance). Only a sufficient reduc- 
tion in the adaptation time r can avoid an instability of 
traffic flow, while large delays in the velocity adjustment 
lead to growing perturbations of traffic flow. 



F Correspondence of the Optimal Velocity 
Model with the Macroscopic Payne Model 

As the Payne model has been claimed to be a macro- 
scopic approximation of the optimal velocity model (see 
Ref. [TU] and citations therein), it is interesting to com- 
pare the instability conditions and characteristic speeds of 
both models. Therefore, let us make the identifications 



P=\> V e (p) = v (- 
d \p 



(111) 
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Then, with the chain rule and the quotient rule of Calculus 
we can derive 



dV e (p) 



dp 



dV c (p) dv Q (l/p) dv Q (d) dd 



dp 
dv (d) 
dd 



dp 



dd 



Inserting this into Eq. ([40]) gives 



dV R 



dp 



1 dv {d) 



dd 



> 



2p c r 



dp 
(112) 



(113) 



or 



dv Q {d c 
dd 



> 



'It 



and 



dV e (p e 



dp 



dVp(de} 

: dd 



.(114) 

where d c = l/p c - This shows the agreement of the insta- 
bility conditions (|40f and (jllpp and of the characteristic 
speeds (|4"Tj) and (|F7)) at the instability threshold. 
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